The advantage of Levy strategies in intermittent search processes 



o 
o 

(N 

Q-r 
D 
C/3 



•4— > 
c3 



i 

C 

O 
o 



> 

(N 

in 
m 

(N 

ON 

O 

o 



X 



Michael A. Lomholt, 1 Tal Koren, 2 Ralf Metzler, 1 and Joseph Klafter 2 

1 Physics Department, Technical University of Munich, D-85747 Garching, Germany 
^School of Chemistry, Tel Aviv University, 69978 Tel Aviv, Israel 

Search strategies based on random walk processes with long-tailed jump length distributions (Levy 
walks) on the one hand and intermittent behavior switching between local search and ballistic 
relocation phases on the other, have been previously shown to be beneficial in stochastic target 
finding problems. We here study a combination of both mechanisms: an intermittent process with 
Levy distributed relocations. We demonstrate how Levy distributed relocations reduce oversampling 
and thus further optimize the intermittent search strategy in the critical situation of rare targets. 
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Random search processes occur in many areas. The 
simplest example is that of passive particles immersed in 
a thermal bath subjecting them to Brownian motion un- 
til encounter, for instance, in chemical reactions [ij. This 
Brownian search dynamics may be accelerated in various 
ways: (i) By a drift toward the reaction center, for in- 
stance, in the time-dependent Onsager problem for diffu- 
sion in an attractive Coulomb potential, or in chemotaxis 
of biological cells 0, 01 • (ii) By combining more than one 
search mechanism due to available interfaces as known 
from gene regulation Q : to find their target sequence on 
the DNA molecule more efficiently, proteins switch be- 
tween 3D bulk diffusion, and ID sliding along the DNA 
(hi) By performing a Levy walk, i.e., a random walk 
whose jump lengths are drawn from a long-tailed distri- 
bution X(x) ~ \x\- 1 - a (0 < a < 2) [f|. The result- 
ing trajectory has diverging variance (a; 2 ) = oo, unless 
a velocity is introduced, and fractal dimension df = a, 
covering space less densely to reduce oversampling, an 
advantage over Brownian search 0, 0, @, Ell EH- 
By intermittent strategics during which local (Brownian) 
search switches with ballistic relocations [12j . 

We here demonstrate for a searcher without orienta- 
tional memory how intermittent and Levy search strate- 
gies can be combined to produce a synergistic strat- 
egy, that for rare targets is more efficient than previ- 
ously introduced intermittent search models. Similarly 
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to Refs. 12 



14 



15j we focus on the ID case, that is 
relevant, for instance, for animals searching for food at 
ecological interfaces (forest edges, coastlines etc.). 

Generalizing the search model from Ref. 14j . we con- 
sider two phases: In phase 1 the searcher looks for the 
target performing diffusive motion with diffusion con- 
stant D. There is a probability per time t^ 1 that the 
searcher leaves this search phase and switches to phase 
2, the relocation phase, where it moves ballistically with 
velocity v in a random direction. The time spent relo- 
cating is drawn from the waiting time distribution ip(t), 
that previously was taken to be exponential (leading to 
a Markovian process) [H , 14 1 , but we relax this assump- 
tion here to show the advantage of Levy strategies. The 
purpose of the relocation phase is to move as quickly as 



possible away from the area that has just been searched, 
and thus the searcher is not scanning for the target in 
this phase. To compare with previous results we take 
a closed cell approach: the search is performed on an 
interval of length L with periodic boundary conditions, 
corresponding to regularly spaced targets with density 
l/L. The model can be formulated as an equation for 
the probability density P(x, t) for the position x of the 
searcher in the search phase: 
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The role of the last term on the right hand side is to 
remove the particle when it arrives at the target placed 
at x = 0. The density pt a (t) thus represents the first 
arrival time at the target, which is determined implicitly 
by the absorbing boundary condition P(x = 0, t) = 0. 
The kernel W(x, t) representing relocations is given by 



W(x,t) 
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^2 S(\x + nL\-vt) 



(2) 



The 5-coupling enforces that the distance traveled in time 
t is vt, and the sum over n renders W(x, t) L-periodic 
in x. ip(t) is related to the spatial distribution of the 
relocations \(x) by ip(t) = 2v\(vt). X(x) is assumed to 
be symmetric around x = (no orientational memory) . 

The search efficiency is quantified by the mean search 
time (t) = J Q dttpfaft). To obtain (t) we Fourier ex- 

pand P(n, t) = J_ L / 2 dx e lknX P(x, t) (n is an integer with 
corresponding wavenumber k n = 2im/L), and Laplace 
transform P(n, u ) = / °° dt e- ut P(n, t), to find 
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uP(n,u) — S n Q — — W(n,u)P(n,u) 

Tl Tl 

1,2 



P(n, u) 

DklP(n,u)- Pla (u). (3) 



The initial distribution is uniform, P(x,t = 0) = l/L, 
since the searcher initially has no information on the po- 
sition of the target. Isolating P(n,u), summing over n 
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1 x 10 Figure 2: Optimal a, and ratio r/ a of search times for optimal 

a versus exponential strategy, as function of L (D — 1 and 
Figure 1: x-t diagram with exponential and Levy relocations, v = 1). All values are calculated using the asymptotic Eqs. 
with n = 37, t 2 = 200, D = 1, v = 0.1 and L = oo. |T0]) and (fl4|) . and corresponding optimal ri and T2. 



(note that X)^L-oo -P( n i u ) = ■P( a; = 0; u ) = 0); an d 
solving for Pf a (u) we find 
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(4) 



In Laplace space the mean search time (t) can be found 
by expanding j»f a at small u since Pf a {u) ~ 1 — (i) it + . . . . 
Be r 2 the average time spent in one relocation event we 
have tp(u) ~ 1 — t 2 u + . . . , and thus arrive at 
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2(n +r 2 ) 



Drxkl + 1 - A(fc„) 



(5) 



Here A(fc„) = W(n,u = 0) = dx e ife " 3; A(a;) is the 
Fourier transform of the relocation length distribution 
at the discrete wavenumbers k n — 2irn/L. We now use 
Eq. ([5]) to determine the search efficiency of (i) Levy and 
(ii) exponentially distributed relocations: 

(i) For Levy distributed relocations we use the sym- 
metric Levy stable law with characteristic function [1 61 ] 



A(fc) 
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2f(l - l/a) 



(6) 



The index a is restricted to 1 < a < 2 so that the mean 
relocation time r 2 is finite. Fig. Q] depicts trajectories for 
cases of exponential and Levy relocations, distinguishing 
the Levy case with its occasional long relocations. 
We introduce three approximations valid for large L: 
(a) Assume that vt 2 S> y/Dr[, i.e., that the mean relo- 
cation distance is much longer than the average distance 
scanned in a typical search phase. We will see that this 
is self-consistent with the obtained optimal values of T\ 
and t 2 that have the same L-scaling for large L. This 
assumption means that Drifc^ and A(fc ra ) are to a good 
approximation non-zero at different n, and we expand 
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D Tl kl + f - A(fc„) Dnh* + 1 l - X(k 



1. (7) 



(b) Assuming that the search range \fDry is much 
smaller than L, we replace the sum over the first term on 
the right hand side of Eq. ((7| by an integral, yielding 
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dn = 



Dnkl + 1 A^fDTl 



(8) 



(c) Approximate the last two terms of Eq. ([7]): as the 
contribution from the singularity at small n dominates 
the sum (note that k n \ n= i — > in the limit of large £), 
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Here C( a ) = S^Li 71 ° 1S ^ ne Riemann C, function. 
Collecting (a) to (c), Eq. §5§ is approximated by 



(t) ~2(r 1 +r 2 ) 
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For honest comparison between Levy and exponential 
strategies, we determine the respective optimal n and r 2 . 
Solving d(t)/dri = and 9(t)/9r 2 = simultaneously, 
we obtain from Eq. (fTU)) that at large L 



Tl 



(6/a a ) 1/(a ~ 1/2) , r 2 ^(6/a^) 1/( ^ 1/2) , (11) 



where (using SI = \/l + 4(a — l)a) 
o=(l + n)/(2[a-l]), 
6 = 2-/D [2a + - 3] C(ar)L a-1 



r(i-a-*) 



(12a) 
(12b) 



such that the optimal t, scale with L like L( a_1 )/( a-1 /2)_ 
According to Eq. (fT0|) . (i) will then scale like 
^(3a-2)/(2a-i)^ implying that for large L the more ef- 
ficient search will occur for a close to 1. However, the 
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Figure 3: Mean search time for Levy (a — 1.35) and exponen- 
tial strategies as function of T2 at asymptotically optimal ri 
(ri = 37.2 for Levy and ri =411 for exponential). We chose 
L = 10 5 , D — 1, v = 1. Simulations versus exact (Eq. ((SJ) 
and asymptotic (Eqs. (1 10 p and (|14[) ) theory. 



prefactor to the L-scaling diverges as a — > 1, so the op- 
timal choice of a will be somewhat larger than 1 for any 
finite L, as demonstrated in Fig. [2] 

(ii) For exponentially distributed relocation with 



(13) 



approximations (a) to (c) also apply, with a — vt 2 . The 
corresponding results for (t) and optimal t% obtain by 
replacing T(l — 1/ot) by 7r/2 and taking a = 2: 
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(14) 
(15) 



These expressions agree with those of Ref . [lj, [17 1 . 

The search time (t) with L for exponential strategies 
scales like L 4 ^ 3 for optimal t\ and r%. This proves that 
Levy strategies with 1 < a < 2 are increasingly more effi- 
cient than the exponential strategies for decreasing target 
density. In Fig. [3] we show (t) as function of relocation 
time t 2 . 

To understand better the a-dependence of the Levy 
strategy we study the first arrival density pi a (t) for large 
L, where again L ^> vt 2 ^> \JDt\. We consider times 
much longer than one relocation-search cycle such that 
ip(u) r-j 1 — t 2 u + . . . , and rewrite Eq. Q as 
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The last expression can be simplified following similar ap- 
proximations as for (t) before. The separation of length 
scales leading to approximation (a) allows us to write 



W (u) 



L 



E 



i 



i 



D Tl kl + 1 nu + 1 - W(n, u) 

(18) 

For the last two terms in Eq. p8[) the contribution at 
small n again dominates the sum (approximation (c)), 
and we expand W(n, u) at small k n and u, finding 
W(n,u) ~ 1 — a a \k n \ a — t 2 u. Collecting the results, 
we see that 



Wb(u)~£ e 
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D Tl kl + 1 (n + r 2 )u + o- Q |fc n |« 

(19) 

We focus on times short enough such that the i- 
periodicity of the problem does not yet play a role, so 
that Laplace space u ^> {cr a \k n \ a \ n=1 ) / {ti + t 2 ). In this 
approximation we replace the sum L _1 Y^=-oo the 
integral dk n /(2Tr), obtaining 



Wb(u) 
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T\l\a sin (rr/a) cr] 
[tt(ri + r 2 )] 1_ 
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(20) 



For shorter times (corresponding to larger it) we discard 
the subdominant second term in Eq. (|20[) . Laplace inver- 
sion of Eq. (fTB|) then produces 



i(n + r 2 ) ' 



(21) 



At later times (smaller u) the second term in Eq. (|2"01) 
dominates, and the plateau (|2ip turns into 



Pfa(t) 
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L(r 1 +T 2 ) 1/a t 1 - 1 /" 



(22) 



The crossover between these two regimes occurs when the 
values of expressions (|22[) and (|2"Tj) become equal, i.e., at 



t~ (n + r 2 ) 



[sm(7r/a)] 2 UT2l a/(a_1) 



(23) 



Note that in Eq. pT]) . 1\JDt\ is the average length 
scanned in a search event. Division by L yields the prob- 
ability to find the target during this phase, and 1/{t\+t 2 ) 
is the rate at which the search phase itself occurs. A 
crucial part in this interpretation is that the probabil- 
ity of searching in a previously scanned area is negli- 
gible. This assumption will break down at some point 
because of the searcher's lack of orientational memory. 
The searcher will then begin to revisit explored regions 
with a reduced probability of finding the target as a re- 
sult. This causes the crossover to the power- law behavior 
Fig. [?] shows the turnover from plateau to inverse 
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— Analytical results (Eq. 21) 




---Analytical results (Eq. 22) 




- Simulation results 











ternatives to Levy strategies, as suggested in Ref. [19j. 
In contrast, the combination of intermittent search and 
Levy relocation strategies turns out to be beneficial. 

Our analysis relies on the assumption that each relo- 
cation is pointed toward a random direction. This will 
be a good model for "non-intelligent" search, similar to 
bacterial movement in absence of chemical or temper- 
ature gradients during which tumbling motion changes 
with directed motion 0]. Intelligent creatures will im- 
prove the target search by partial or complete memory, 
avoiding previously visited locations. It will be interest- 
ing to study in more detail models with search memory. 
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Figure 4: First arrival density versus time. The crosses are 
simulation data, while the straight lines are the intermediate 
regimes of Eq. (|21[) and Eq. (|22|l . Parameters are: ri = 35, 
r 2 — 50, L — 10 4 , a = 1.75, v = 1 and D = 1. 

power-law of the first arrival. At even longer times, finite 
size effects cause a turnover to an exponential decay. 

From Eq. (ffi?)) the advantage of having a close to unity 
at large L becomes evident: the presence of rare but long 
relocation events reduces the risk of rescanning already 
visited areas which will be important for large L. How- 
ever, the downside to choosing an a- value too close to 1 
is that an increased amount of very long relocations im- 
plies an increased amount of very short ones too, as the 
average distance is fixed by vti [l8[. This means that 
the crossover to the less favorable situation described by 
Eq. (f!?2)) happens earlier, so that larger a becomes more 
efficient for shorter search times relevant at smaller L. 

From a more general perspective, intermittent strate- 
gies are beneficial when purely diffusive search would slow 
down over time due to the increasing number of returns 
to previously scanned areas (oversampling) . Choosing an 
exponential strategy for the relocation events, however, 
only partially solves this problem, as this strategy is still 
governed by the central limit theorem (CLT). Thus, the 
problem of oversampling merely becomes postponed to 
later times. Conversely, Levy-intermittent strategies are 
not bound to the CLT, rendering them advantageous in 
the search for rare targets. Although less pronounced, 
the problem of oversampling still occurs in two dimen- 
sional search studied in (l9j ]. Thus, Levy strategies are 
expected to be advantageous in this well. 

We have shown that for intermittent search strate- 
gies Levy distributed relocations are advantageous over 
exponential distributions when targets are sparse, be- 
cause rare long relocations reduce the eventually occur- 
ring problem of oversampling. Thus we advocate that 
intermittent strategies should not be thought of as al- 
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